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The study of the entanglement entropy and entanglement spectrum has proven to be very fruitful 
in identifying topological phases of matter. Typically, one performs numerical studies of finite- 
size systems. However, there are few rigorous results in this regard. We revisit the problem of 
determining the rank of the “particle entanglement spectrum” of the Laughlin states. We reformulate 
the problem into a problem concerning the ideal of symmetric polynomials that vanish under the 
formation of several clusters of particles. We introduce an explicit generating set of this ideal, and 
we prove that polynomials in this ideal have a total degree that is bounded from below. We discuss 
the difficulty in proving the same bound on the degree of any of the variables, which is necessary to 
determine the rank of the particle entanglement spectrum. 


I. INTRODUCTION 

The study of topological phases of matter has benefited greatly from consideringthe entanglement 
properties of the ground states of topological phases. The work of Kitaev and Preskill 1 ! and of Levin and 
WerP revealed that the entanglement entropy is a good probe of the topological nature of a system and 
provides a measure for the particle content of the topological phas^.The entanglement entropy of a pure 
quantum state |T) relative to a bipartite partition of the total Hilbert space TL = TLa < 8 > 'Hb provides a 
measure of the entanglement of |T). The entanglement entropy is defined as the Von Neumann entropy 
S of the reduced density matrix of either one of the two parts. For instance 

S =-Tr(p A log p A ), (1) 


where p A = Tr# (|\I/) (\E'|). 

In the context of the fractional quantum Hall (FQH) effect, various ways to partition the Hilbert space 
were proposecP. Of particular importance, is the spatial partitioning scheme in which the system is split 
into two regions A and B separated by a real-space cut of length L. For a system exhibiting topological 
order the real-space entanglement entropy is of the form 2 3 

S = aL - 7 H-, (2) 

where • • • stands for subdominant terms as L becomes large. The subdominant term 7 is universal, 
and depends only on the nature of the topological phase. It bears the name topological entanglement 
entropy , and is a measure for the particle content of the topological phas^. The first term aL, while 
non-universal, means that the amount of entanglement is proportional to the length of the boundary 
separating the two regions. This property called area law has appeared in various areas of physics, 
such as black-hole physics and quantum information. For a quantum many-body state, this property is 
of particular importance since it opens the way to extremely efficient numerical simulations such as the 
Density Matrix Renomalization GroupPand Matrix Product Stated methods. For FQH state this avenue 
of research was successfully undertaker^®! anc j opened the way to a reliable microscopic calculation of 
quasi-holes properties such as radius and braiding^ 0 . 

Although the real-space cut is of paramount importance in the study of topological phases of matter, 
there are other natural ways to partition a quantum Hall system: the orbital cut, and the particle cwPl 
While, in principle, the entanglement entropy behaves according to the area law Eq. |2]) only for real- 
space cuts, it was numerically observecP 1 that the area law is also valid for orbital cuts. In this paper 
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we will concentrate on the particle cut, in which one numbers the (identical) particles constituting the 
phase (for instance, the electrons in the quantum Hall case), and one declares the particles numbered 
1,2,..., Na to belong to subsystem A, while the remaining particles numbered Na + 1, Na + 2, ..., 
N belong to subsystem B. The spectrum of the reduced density matrix obtained by tracing out the 
particles in subsystem B is the “particle entanglement spectrum” (PESj-3 

While the entanglement entropy S provides a good probe of topological order, the topological entan¬ 
glement entropy 7 does not determine unambiguously the universality class of the topological state. Li 
and Haldane 13 realized that the spectrum of Ha = — log pa itself contains much more information than 
the entanglement entropy. They proposed to use the low lying part of this entanglement spectrum as a 
“fingerprint” of the topological phase. To be more specific, under a bipartition H = Ha < 8 > H B , a pure 
quantum state |\f r ) admits a Schmidt decomposition 

w = E e_c</2 i^>®i^)> ( 3 ) 


where the e ^/ 2 ’s are positive numbers called the Schmidt singular values, while \ipf) and \ipf) form 
orthonormal sets in Ha and H B , respectively. The reduced density matrix is then simply 

PA = E e_?i IV’i l }(V , i 4 | =e~ HA . (4) 

i 

The entanglement spectrum is the set of all entanglement energies The bipartition can be chosen 
to preserve as much symmetry as possible, which in turn yields quantum numbers for the £j’s, such 
as the momentum along the cut. Li and Haldane observed that-per momentum sector-the number of 
entanglement energies reproduces exactly the number of gapless edge modes. They proposed that tracing 
out the degrees of freedom of part B introduces a virtual edge for part A. The Li-Haldane conjecture is 
therefore two-fold. For a FQH state in the thermodynamic limit: 

(i) the entanglement energies and edge modes have the same counting, 

(ii) the entanglement spectrum is proportional to the (edge) CFT spectrum. 

It is now understood that (ii) can only hold in the case of a real-space cut, which maintains locality 
along the cufP^ESl. p or an orbital cut the entanglement Hamiltonian Ha has no reason to be local. On 
the other hand the point (i) holds irrespective of the cut for model wave functions that can be written as 
correlation functions in a CFT. Such wave functions are precisely of MPS fornP^, and the CFT Hilbert 
space provides a one-to-one mapping^ between edge modes and entanglement energies. 

While the agreement between the counting of the number of modes in the entanglement spectrum 
and the counting of the edge excitations is well understood, in practice, this fingerprint is used for finite 
-size systems. The entanglement counting develops finite-size effects, which naively have no structure. 
However, it has been conjectured and numerically substantiated ^ 9 that there is a counting principle 
underlying the finite-size entanglement counting of model states. Before we turn to the main topic of 
this paper, namely the PES for the Laughlin states, we mention that the entanglement entropy and 
spectra are active studied, see references^^for some studies in the context of the quantum Hall effect. 

We now consider the ground state T) of a model quantum Hall state, such as the Laughlin 27 or 
Moore-ReacpSl state, that are the exact zero energy states of a model Hamiltonian. For a given number 
N of particles, this is the unique zero-energy state of a model Hamiltonian that occurs at the following 
number of flux quanta 


Nt = ^N~S, (5) 

where v is the filling fraction and S is the shift. Now suppose that the N particles are divided into 
two groups, group A containing Na of the particles, and group B containing the rest N B = N — Na 
of them. Without any loss of generality, we can assume that Na < N B . Let x = (aq, X 2 , • ■ • , Xn a ) and 
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V = (z/iiZ/ 2 !'' • ,Vn b ) be the coordinates of particles in A and B, respectively. The Schmidt decomposi¬ 
tion of the wave function '!'( 2 ), that is 

^( z ) = J2 e ~ ii/2 ' l Pi L ( x )' l l’i 3 (y)’ ( 6 ) 


involves wave functions i/>f(x ) for Na particles. After tracing out the particles of part B, we are left 
with a reduced system of Na particles, but the amount of flux remains the same, namely 

JV* = ^N a - S + AN+, (7) 

where A = v~ 1 Nb■ The presence of this excess flux ATV^ indicates that we should view the reduced 
system as one with Na particles, in the presence of quasi-hole excitations. For a real system with Na 
particles, and Aexcess flux quanta the number of zero-energy states of the model Hamiltonian (which 
we will call the number of quasi-hole states) can often be obtained exactljS2H23) p or instance, in the case 
of the v = 1/m Laughlin case, quasi-hole states of Na particles in A N^ excess flux are of the form 

= P z (x)Y]_(xi - Xj) m , (8) 

i<j 

where Pi(x) is a symmetric polynomial in Na variables with degree in each variable at most A N^. The 
number of quasi-hole states is therefore 


This number forms an upper bound for the rank of the reduced density matrbPI - 11 23 ^. 

From numerical investigations, it is known that in all cases considered, this upper bound is in fact 
reached^). This observation has led to the “rank saturation” conjecture, which can be thought of as a 
finite-size version of the Li-Haldane conjecture, namely, the entanglement level counting of the PES of a 
model state is equal to the number of bulk quasi-hole states. This means that the states ipf(x) appearing 
in the Schmidt decomposition of 'S’(z) span all the quasi-hole states of Na particles in AN^ excess flux. 
Proving analytically that this upper bound is indeed reached has proven to be a difficult problem. 

In this paper, we revisit this problem for the general v = — Laughlin states. We start by considering 
the v = 1 Laughlin state, which is simply the Slater determinant of the completely filled lowest Landau 
level. We explain how to obtain the rank of the reduced density matrix of the particle entanglement 
spectrum in this case. To do so, we will make some use of the properties of symmetric polynomials. 
To get a grip on the v = A Laughlin states, we then use the following strategy. After partitioning the 
particles into two sets A and 5, we “split” the Nb particles in part B into mNs particles, and consider 
the v = 1 Laughlin state of the system thus obtained. For this system, we already obtained the rank 
of the reduced density matrix. If one can show that clustering the mNs particles into groups of size 
to, does not lead to a smaller rank of the reduced density matrix, one deduces the rank of the reduced 
density matrix for the v = -b Laughlin state, and shows that the upper bound is indeed reached. 

The hard step in the strategy outlined above is to show that the clustering of the mN B particles 
into Nb groups of of size to does not reduce the rank of the reduced density matrix. Proving this 
statement turns out to be highly non-trivial. As we explain in the main text, one has to show that 
there is no (non-zero) symmetric polynomial in mNs variables that vanishes under the formation of Nb 
groups of variables each of size to, and whose degree in any of the variables is Nb or less. Although 
we did not fully succeed in proving this statement, we did make substantial progress. In particular, we 
constructed an explicit generating set of the ideal of polynomials that vanish under this clustering. Using 
this construction we were able to show that a non-zero symmetric polynomial in mNs variables that 
vanishes under the formation of Ng groups of variables each of size to must have a total degree at least 
Nb + 1. Proving this weaker statement is already a non-trivial result, mainly because the positions of 
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the various clusters can be arbitrary, which means that the clustering condition is non-local. Moreover 
we were able to prove that all polynomials in the generating set have a degree at least Nb + 1. 

The outline of the article is as follows. In section [TlJ we introduce the notion of partitions, and 
several types of symmetric polynomials, that we make use of throughout the article. The PES of the 
v = 1 Laughlin state is discussed in section |III| We continue in section [IV] by explaining how the result 
for v = 1 can be used to make progress for the v = 1/m Laughlin states, and recast the problem in 
terms of clustering properties of symmetric polynomials. In section [Vj we prove that the total degree of 
polynomials that vanish under clustering is bounded from below, and provide an explicit construction for 
such polynomials in general. In section |VII[ we make some comments on why it is much harder to prove 
that not only the total degree, but also the degree of any variable for polynomials that vanish under the 
clustering is bounded from below. In addition, we provide a proof for the statement in the case where 
one forms two clusters of size m. Finally, we discuss our results in section [VIII| In the Appendix A, we 
derive some properties of the polynomials which are used in section |VJ and in Appendix B we provide 
an alternate set of polynomials that can be used in the proof of section [V| 


II. SOME NOTATION 

In this section we introduce some definitions and notations that are used in what follows. We start by 
introducing the notion of partitions, which play a central role in the theory of symmetric polynomials. 
For a general introduction to the subject of partitions, we refer tcP^ and for the theory of symmetric 
polynomials tcP®. 


A. Partitions 

For a positive integer V, a non-increasing sequence A = (Ai, A 2 ,..., A r ) of strictly positive integers Ai, 
A 2 , ..., A r is called an r-partition of V if A,; = V. The A/s are the parts of A, and r is called the 
length of A, which is denoted by l( A). We call V the weight of A, which is denoted by |A|. We write 
A b V to indicate that A is a partition of V. By convention, A = 0 is the only partition of zero which 
we call the empty partition. The number of parts of partition A which are equal to a given integer j is 
denoted by rij{ A) or simply rij. We also define 
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( io ) 

j=i 

Finally, The set of all partitions of V is denoted by Par(T>). It is not too hard to convince oneself (see 
Ref. (35 J that the number of partitions with at most r parts and each part at most d is equal to ( r + d ). 


B. Symmetric polynomials 

In what comes, we will be dealing with the ring Atv of symmetric polynomials in N variables. A 
polynomial P is called a symmetric polynomial in N variables, if for all permutations a of {1,..., N}, 

P = P(X 1,...,X N ). (11) 

The degree d of a symmetric polynomial is simply the degree in one of its variables. 

A polynomial P(aq,..., Xjy) is called homogeneous of total degree V , if for any real number l, 

P(l X 15 • • • , l Xjv) l P{x 1 , . . . ,^tv)- 


( 12 ) 
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For instance, the polynomial P(x i,x 2 ) = x\x 2 +X\x\ is a homogeneous symmetric polynomial of degree 
cL = 2 and total degree V = 3. 

There are different bases that one can consider for An. A natural one, is given by the set of, so-called, 
symmetric monomials. Given a partition A = (Ai,A 2 ,--- , A r ) with r < N, the symmetric monomial 
m\(x\, ■ ■ ■ , Xn) is defined as 

m\(x,X N ) ; = 5Z*a(l)*a(2) ••• X a( r ) a: °(r+l)"' a; °W’ ( 13 ) 

a 

where the sum is over all distinct permutations a of the parts of A, and it is defined to be 1 if A is the 
empty partition. On the other hand if r > N we set m\{x 1 , • • • , Xn) = 0. For example, 

m( 2 ,i,i)(xi,x 2 ,x 3 ) = x\x 2 x 3 + xix\x 3 + x\x 2 x\, 
m( 2 , 1 ) (xi, x 2 ,x 3 ) = x\x 2 + x\x 2 + x\x 3 + x\x\ + x\x 3 + x 2 xl, (14) 

ixi( 2 ,\,\){x\,x 2 ) = 0 . 


When studying rank saturation of the PES for the Laughlin state, finite-size effects imply an upper 
bound for the degree of polynomials. We will therefore be led to consider the space of symmetric 
polynomials in N variables, with degree (in each of the variables) at most d. A basis for this space is 
given by the symmetric monomials m\(x 1 ,..., Xn) corresponding to partitions A with at most N parts 
and each part at most d. Therefore, we have 

dim (A d N ) = (^ + d ' S j. (15) 

Another important family of symmetric polynomials is the set of elementary symmetric polynomials. 
The elementary symmetric polynomials that are labelled by an integer n are defined in terms of symmetric 
monomials as e n := rn ( i 2 )- For instance, 

eo(xi,x 2 ,x 3 ) = 1 , 
e 1 (x 1 ,x 2 , x 3 ) = X! + x 2 + x 3 , 
e 2 (xi,x 2 , x 3 ) = xix 2 + 2 : 1 X 3 + x 2 x 3 , 
e 3 (xi,x 2 ,x 3 ) = xix 2 x 3 , 

e n > 4 (xi,x 2 ,x 3 ) = 0 . (16) 

For a partition A = (Ai,..., A r ), the elementary symmetric polynomial e\ is defined as e\ := • • • e\ r . 

As an example, 


e ( 2 ,i,i)(xi, x 2 ) = e 2 (x!, x 2 )e 1 (x 1 ,x 2 )ei(x 1 , x 2 ) 

= X\X 2 (X\ + x 2 ) 2 . 

It is known that the set of all polynomials e\(x\,... ,xjv), where A is a partition with at most d parts 
and each part at most N, forms a basis of the space A^. 

Lastly, the power sum symmetric polynomials , defined as 


Pi(xi, ... ,x N ) := x\ H-b x)v , (18) 

are of special importance. In fact, the set {pi,p 2 ,... ,Pn} generates Ajv- This means that any symmetric 
polynomial P in N variables can be written as a polynomial in (pi,... ,pn)- In other words, the set of 
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polynomials p\ := p\ x ■ ■ ■ p\ r , where A is a partition with each part at most N, forms a basis of Ajy. For 
example, 

e 2 = XiX i = Pl o ^ > ( 19 ) 

i<j 

independent of the number of variables N. Most importantly, the decomposition of any symmetric 
polynomial P in N variables as a polynomial in (pi,...,pjv) is unique. One should note that, unlike 
for the symmetric monomials m\ and the elementary symmetric polynomials e\, there is no natural 
restriction on A such that the corresponding p\ s form a basis for Afj. 


III. THE v = 1 STATE 


To obtain the rank of the reduced density matrix of the Laughlin states in the case of the “particle 
cut”, we start by considering the simplest case, the v = 1 Laughlin state, which is just a single Slater 
determinant, 

'IwO z N )= (Zi-Zj), (20) 

up to a geometry-dependent Gaussian factor. For instance, the plane and sphere geometry give rise to 
different Gaussian factors, inherited from the respective inner products. However, for our purposes the 
precise form of the Gaussian factor is irrelevant. The results presented in this paper involves only the 
notion of linear independence, and does not refer to the notion of orthogonality. As a consequence, the 
underlying inner product plays no role and our result is equally valid on the plane, sphere, and cylinder. 

Now suppose that the N particles are divided into two groups A, containing Na of the particles, and 
B containing Ng = N — Na particles. At this stage we do not assume Na < Ng ■ Let us rename the 
coordinates of particles in A and B to x = (xi,X 2 , ■ ■ ■ , Xn a ) and y = (jq, j/ 2 , • ■ ■ , VNb)i respectively. The 
rank of the reduced density matrix in the case of such a particle cut can be obtained from a decomposition 
of the wave function ^(z) of the form 




0 = 


( 21 ) 


where the set of wave functions (resp. ijjf) are independent. Note that this is not quite a Schmidt 
decomposition since we do not demand the s to form an orthonormal set. Although this is not a 
Schmidt decomposition, the number of terms in the sum is equal to the Schmidt rank, or equivalently, 
to the rank of the reduced density matrix. Therefore, we will call the decomposition (211 a Schmidt 
decomposition, although strictly speaking this is an abusive notation. 

Before we explicitly write the v = 1 Laughlin state in such a “Schmidt-decomposed” form, we note 
that we can obtain the rank of the reduced density matrix in the v = 1 case in a straightforward way. 
This state is simply obtained by filling the Landau orbitals from 0 up to A$ = N — 1 

|^=i) = |111 - • • 111). (22) 

The Schmidt decomposition relative to particle cut amounts to choose Na out of the N particles 


|^=i) oc I 111 • • • 1100 • • • 0) <g) 1000 • • • 0011 • • • 1) 

V -v- y 

N a 


+I 111---10 10-- 

Na 


|000 • • • 0101 - - - 1 ) + ■■• 

+ I 000- -00 11 • • • 1) <g) |111 • • • 1100 • • • 0), (23) 

N b 
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which means that the rank of the reduced density matrix is given by the number of ways in which the 
Na particles of system A can be divided over the number of orbitals. The number of orbitals is given 
by N = Na + Nb, so we obtain that the rank of the reduced density matrix is given by ( Na ^ b )- We 
note that the same result can be obtained directly from the wave function 1 ^ which is a single Slater 
determinant T 1/= i (z) = (N ~ z i)- 

It is instructive to perform a more explicit Schmidt decomposition of the v = 1 Laughlin wave function. 
We start by writing the state explicitly in terms of the variables Xi and yi of groups A and B , respectively. 
Dropping the exponential factors, we have 


^u=i (z) = ^ v =i{x) 



-yj) 'iw(y)- 


(24) 


We are going to use the following result!^ 


N a N b 

n n^ 1 + ^yi)=n m\{x)e\{y). (25) 

i=1J=1 A 

where the sum is over all partitions A with maximally Na parts, and each part being maximally Ng, 
i.e., all partitions which fit in a rectangle of height Na and width Ng- Thus, 

_ Vi) = H(“ 1 ) |A ' m A( a; ) e A(y)- ( 26 ) 

i,j A 


Here, we used the relation 


IK* m A (-l/a;) 




(27) 


where the partition A stands for the complement of A with respect to the rectangle of height Na and 
width Ng- In addition, (—1/a:) is shorthand for (— l/x±, ..., 1/xn a ). As an example, it is shown in 
Fig. 0 that for Na = 3, Ng = 4 and A = (2,1), one finds A = (4,3,2). We then obtain a Schmidt 


N B = 4 


A = 

(2,1) 





A = 

(4,3,2) 






FIG. 1: The relation between the partition A and its complement A for given Na and Nb- 
decomposition for the v = 1 Laughlin state 

^ =1 (z) = ^(-l)' A '^(a:) g f(y), 

A 


(28) 
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where 


Qx(x) = mx(x)V v= i{x), 


(29) 


and 


Qxiv) = e\(y) ^v=i{y)- 


(30) 


A few remarks on this formula are in order here. The number of terms in the sum is most important 
here. The sum over A, is over all partitions with maximally Na parts, and each part being maximally 
Nb, i.e., all partitions which fit in a rectangle of height Na and width N b . There are ( Na ^^ b ) such 
partitions. The rank of the reduced density matrix of the v = 1 Laughlin state is thus given by 


fN A + N B \ 

V N A )' 


(31) 


and we recover the dimension of . It is straightforward to check that this is also the dimension of 
the set of anti-symmetric polynomials in Na variables with maximum degree N = Na + N b , which is 
nothing but the space of ‘quasi-hole’ states for the non-interacting v = 1 case. The set of polynomials 
q ^ (resp. qB) forms a basis for the space of anti-symmetric polynomials in Na (resp. N b ) variables with 
maximum degree N. Note that this result is symmetric under exchange of A and B , and in particular 
it holds whether or not Na < N b . This is a particularity of the v = 1 case and it will no longer be true 
for the v = N Laughlin state with m > 1. 


IV. SCHMIDT DECOMPOSITION OF THE v = A LAUGHLIN STATE 

m 

We are now going to compute the rank of the reduced density matrix for the generic v = 1/m Laughlin 
state 


(Zi - Zj) m . 

l<i<j<N 


(32) 


As usual we divide the particles into two groups A and B, containing Na and N b = N — Na of them, 
respectively, and we assume that Na < N b . We are interested in obtaining a Schmidt decomposition of 
this state. As for the v = 1 case, we can write 


Na N b 

Vm(z) = 'km(x) ( IJ(aq - yj) m | 

,*= 1.7 = 1 


(33) 


Proving that the rank of the PES for the v = 1/m Laughlin state is saturated boils down to finding a 
Schmidt decomposition for the wave function ( |33[ ) and counting the number of terms in the decomposition. 
As in the v = 1 case, one need to take care of only the middle term of the wave function (331 


N A Nb 

: =nil 

*= 11=1 


(34) 


To do so we start from 


Na mN B 

$i (x,w )=n n ( Xi - 

i=i j=i 


(35) 
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for w = (w 1 ,.. .,WmN B )- From (26), 


$i(*,w) = ^(-l) |A| mx(a;)eA(w), 


(36) 


but this time the sum is over all partitions A which fit in a rectangle of height Na and width mN B . 
We then relate <f> m to •!>! through the clustering transformation, which is a linear transformation from 
A m Ar B to An b defined as follow. 

To a symmetric polynomial P(w) in mN B variables, we associate the polynomial variables C m (P(y)) 
in Nb variables obtained by regrouping the particles into clusters of to, i.e., 


(C m P)(yi, 2 / 2 , • ■ •, vn b ) = P(yi,yi, ■ ■ ■, yi, 2 / 2 , 2 / 2 , • • •, 2 / 2 , ■ • •, vn b , vn b ,...,y NB )- 


(37) 


It is easy to see that after clustering <f>i becomes 3> m , i.e., 

Cm 1 ) = ^ m • 


(38) 


Applying the clustering transformation to both sides of Eq. (361 results in 


®m(x,y) = ^2(-l) lxl mx(x)C m (ex){y)- 


(39) 


As mentioned earlier, the sum is over all partitions A with maximally Na parts, and each part being 
maximally mN B ■ There are { Na ~n™ N b ) such partitions. This is precisely the number of Laughlin quasi¬ 
hole states for Na particles in AN$ = mN B extra fluxes, and we recover the usual upper bound for the 
rank of the reduced density matrix. 

Rank saturation of the PES for the v = 1 /to Laughlin state boils down to the following, non-trivial 
result: the polynomials C m (ex) are independent. More precisely, one has to prove that the linear trans¬ 
formation 


A Na 

Iy mN B 


A mN a 
1 V N b 


P{w) - > Cm{P){y) 


is injective as long as Na < N B . Since dim N^ A > dim A A , it is sufficient to show that this linear 
map has a trivial kernel. Namely, besides P = 0, no polynomial in mN B variables and maximum degree 
Na can vanish under the clustering transformation. 


V. CLUSTERING PROPERTIES OF SYMMETRIC POLYNOMIALS 


In this section we are going to describe the ideal of symmetric polynomials in q = m N variables that 
vanishes under the clustering transformation (37). In particular, we are going to construct a generating 


set of this subspace, and prove that a non-zero symmetric polynomial in q = to N variables that vanishes 
under the clustering transformation has a total degree V of at least N + 1. We are also going to prove 
that this symmetric polynomial of minimal total degree is unique (up to a scaling numerical factor). 

The statement that the total degree of a symmetric polynomial in q = mN variables that vanish 
under the clustering conditions is at least N + 1 is a weaker statement than stating that the degree of 
each variable is at least N + 1, but easier to prove. After finishing the proof of the statement on the 
total degree, we comment on how one might prove the stronger statement, limiting the degree of the 
polynomials. 

As a warmup, we start with two simple examples, which we will come back to after the proof. We 
start with the case m = 2 and N = 1, i.e., we are looking for a symmetric polynomial in two variables 
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yi and 2/2 , of total degree 2, that vanishes when yi = y 2 . It is easy to see that the polynomial has degree 
two, namely (y 1 - y 2 ) 2 ■ 

The case m = 2 and N = 2 is already more complicated. With some thought, one can construct a 
total degree 3 symmetric polynomial in four variables y-\,..., j /4 , that vanishes when yi = y 2 and y 3 = j/ 4 , 
namely ( yi + y 2 - y 3 ~ Vi) ( 2/1 - 2/2 + 2/3 - 2/4) ( 2/1 - 2/2 - 2/3 + 2 / 4 )- It is already less trivial to convince 
oneself that no lower degree symmetric polynomial with the same vanishing properties exists. Upon 
increasing the m and N, even finding polynomials with the correct vanishing properties becomes a hard 
problem, which is caused by the non-locality of their defining property. Namely, polynomials have to 
vanish, independent of the position of the various clusters. As we indicated above, we solve this problem 
in a constructive way. 

Our construction is motivated by the following observation. The ring A q of symmetric polynomials in 
q = mN variables is generated by {pi,... ,p q }, and the power sum polynomials Pi have a very simple 
behavior under the clustering transformation ( |37| ), namely 

Cm{Pi) =mpi. ( 40 ) 


However, after the clustering transformation there are only N variables left. This means that 
An = {pi, ■ ■ ■ ,Pn} forms a minimal set of generators, and the polynomials {pn+ 1 , ■ ■ ■ ,p q } are no longer 
independent after being clustered. The generators pi are not very convenient to describe the clustering 
transformation, and this is why we introduce a new set of generators A q = { p \,... ,p q } of A q as 


Pn = 



(41) 


Alternatively, the polynomials p n can be defined in terms of the polynomials rof Appendix |a| through 
The main property of these new polynomials is that they behave nicely under clus- 


„ _ (zUTJ 1 /™ 

Fn — n\ n 

tering: 


Cm(pn) = e„, n = l,...,AT, (42) 

Cmipn) = 0 , n>N , (43) 


as inherited from the properties of rn '* described in Appendix [a| We also introduce the p\ in the usual 
way as p\ = p\ 1 .. .p\ r , where r is the length of A. In terms of these modified power sums p n , it is now 
relatively easy to describe the ideal of polynomials in A q that vanish under the clustering transformation 

Theorem 1. The set {^1^+1 < Ai < q} forms a basis of the ideal of symmetric polynomials in q = mN 
variables that vanish under the clustering transformation. In other words, this ideal is generated by the 
set {pn+i,Pn+ 2, ■ ■ ■ ,Pq}- 


Proof. Suppose that P is a symmetric polynomial in q variables. Because A q is a generating set, there 
exists a polynomial R in q variables such that 

P = R(Pl,---,Pq)- 

Generically, there are two kinds of monomials in the polynomial R. Those that depend only on the first 
N variables Pi,... ,Pn, and the ones that depend on at least one of the p n , with n > N. Accordingly, R 
can be decomposed uniquely into a sum of two polynomials 

R(,Pi, ■ ■ ■ ,Pq) = A{p- l, ... ,p N ) + B(pi,... ,p q ). 

Thus, by construction, C m (B ) = 0. It is now straightforward to check that C m (i?) = 0 if and only if 
A = 0, since 


C m (P) = A(e e N ), 
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and the {ei,..., ex} are algebraically independent in N variables. Therefore the set {pa| 77+1 < Ai < q} 
is a basis of the kernel of the clustering transformation. ■ 

Corollary 1. The only symmetric polynomial P in mN variables with total degree N or less that vanishes 
under the clustering transformation is P = 0. Moreover, pn+i is the unique (up to an overall factor) 
symmetric polynomial in q variables and total degree N + 1 that vanishes under this clustering. 

Since the modified power sum p n has total degree n, this corollary follows directly from Theorem 1. 
Let us illustrate this with two simple examples. 

Example 1. Consider the simplest non-trivial case where q = 2, N = 1, and m = 2. In this case, the 
clustering condition is y\ = 2/2 • The definition of p 2 yields 

P2 = ^ (pi ~ 2 p 2 ) = ~\(yi ~ v?) 2 ’ ( 44 ) 

which reproduces the expected result. 

Example 2. As another example, let q = 4, N = 2, and m = 2. This time, clustering conditions are 
yi = y 2 and 2/3 = j/ 4 . We have 

P3 = ^{pl ~ &V 2 P 1 + 8 p 3 ). (45) 

For q = 4 variables, this is 

P3 = yx(2/i + 2/2 - 2/3 - 2 / 4 )( 2/1 - 2/2 + 2/3 - 2 / 4 )( 2/1 - 2/2 - 2/3 + 2 / 4 )- (46) 

16 

We should note that these two examples are not representative for the general case, in the sense that 
the polynomials Pn+ 1 do not generically factorize to a simpler form. For instance, for N = 3 and m = 2, 
we have 

P4 = - 12 P2P\ + I2pl + 32pip 3 - 48p 4 ), 

which does not have a simple factorized form when restricting to q = 6 variables. 

Conjecture 1. There is no non-zero symmetric polynomial P in mN variables with degree N or less 
that vanishes under the clustering transformation. 

While we know that the modified power sum p n has degree n (see Appendix |A|, this is not sufficient 
to prove this conjecture. 


VI. 527(2) INVARIANCE 


In the context of the fractional quantum Hall effect, there is a natural action of 517(2) on A^7 coming 
from the rotational invariance of the sphere. The angular momentum operators on the sphere 36 are 


zr 

L 3 


L+ 



(47) 

(48) 

(49) 
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Every polynomial P in A 1 ^ has a SU(2) symmetric f2(P) with opposite angular momentum L 3 given 

by 


n(P)(zi,--- ,Z N ) 





(50) 


Under this Z 2 operation, L and L + are exchanged and L 3 —► —A 3 . 

These linear operators are compatible with the clustering, in the sense that C m L l = L l C m . Note that 
in these identities the SU(2) operators in the l.h.s. act in A^jy, while in the r.h.s. they act in Aff NiS : 


A N * 
ly mN 


, A N® 
^ ly mN 



a mNq> 
1V N 


a™ n * 


(51) 


The same is true for the Z 2 operation Q. These commuting properties are straightforward to check for 
A 3 , L~ , and for ft. Therefore, it also holds for L + = flA _ IA For instance, the clustering transformation 
clearly preserves the total degree, hence the action of clustering commutes with L 3 . Likewise, A” being 
the generator of global translations, it commutes with the clustering. The following theorem follows 
immediately: 

Theorem 2. The ideal of symmetric polynomials in mN variables that vanish under the clustering 
transformation is invariant under the action of L l and f2. 

Corollary 2. The polynomial pn +i is translationally invariant. 

The polynomial Pn+i vanishes under clustering, and therefore so does L~pjy + i. If this last polynomial 
was non-zero, it would have a total degree N, which is forbidden by Corollary 1. Therefore L~pn+i = 0 
and pn+i is translationally invariant. 

In fact, it is possible to directly calculate L~p i: and one gets 

L~pi = (N + 1 - i)pi-i. (52) 


(x) 

Note that this results only hold for q = mN variables. This follows from the behavior of r n under 
translations, which is given in Appendix |A 4| 

Since the kernel of C m is invariant under the action of L 1 , it can be decomposed into irreducible 
representations of SU(2). In order to prove that non zero polynomials that vanish under the clustering 
have degree at least N + 1, it is therefore sufficient to prove it for lowest weights, that is to say translation 
invariant polynomials. Therefore Conjecture [l] is equivalent to the following: 

Conjecture 2. The only translationally invariant symmetric polynomial P in mN variables with degree 
N or less that vanishes under the clustering transformation is P = 0. 


VII. A POSSIBLE ROAD TOWARDS FINISHING THE PROOF 


As we saw in the previous section, we were able to prove that the total degree of a symmetric polynomial 
is at least IV+ 1, if the polynomial vanishes under the clustering transformation (371. However, we would 
like to show that the the maximum degree of any of the variables (i.e., the number of fluxes IV^) is at 
least N + 1. Proving this statement turns out to be much harder than it looks at first. One of the reasons 
is that the clustering we consider is a non-local process. Namely, the positions of the various clusters 
are arbitrary. Therefore, proving that the total degree is bounded from below is already a non-trivial 
result. What makes proving a bound on the total degree more tractable in comparison to proving a 
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bound on the degree, is that upon taking linear combinations, the total degree of the polynomials does 
not change, provided the resulting polynomial does not vanish. The degree of the polynomial, however, 
can be lowered by taking linear combinations. 

To show that the rank of the reduced density matrix for the particle cut does indeed satisfy the upper 
bound given in the introduction, it suffices to prove that the clustering map C m : A^jy — ► A’^' 4 , is 
injective if Na < Nb- In the case Na = Nb, the map C m would then actually be bijective. One possible 
route in trying to prove this, is to find two suitable bases for A^ N and A™ N in which the map C m acts in 
an upper-triangular way, and then check that all diagonal elements are non-zero. We did not, however, 
succeed in finding suitable bases. 

A completely different route to prove that the rank of the reduced density matrix is given by the 
upper bound is to try to make use of the results for the Read-Rezayi states 3 ^. These states are defined 
by the property that they vanish if k +1 particles are put at the same location (in their simplest bosonic 
incarnation). In particular, it is known exactly that how many symmetric polynomials there are, that 
satisfy this clustering condition, for an arbitrary number of particles, and arbitrary degre^- 32 33 4 In 
addition, there are explicit expressions for these polynomial^! see alstPS Using these results, we can 
prove the wanted result for N = 2 and arbitrary to. That is, we can show that any symmetric polynomial 
in 2m variables, that vanishes if two clusters of to variables each are formed, has degree at least three. 

To do so, assume that P is a polynomial in 2 to variables that vanishes under the clustering, C m P = 0. 
We know that the total degree of this polynomial is at least 3, and we want to show that the minimal 
value of the degree is three as well. 

To show this, we note that the polynomial P also vanishes if we make one big cluster of 2 to variables. 
From the results on the Read-Rezayi states, we know that such symmetric polynomials have degree at 
least two (it vanishes, so it should vanish quadratically), and that it is unique (up to an overall factor). 
In addition, we know an explicit form of this polynomial P ', namely 


P'{zi, Z 2 m) = S[{zi - z 3 )(z 2 - z 3 )], (53) 

where S denotes the complete symmetrization over all 2 to. variables. In this case, by inspection, one can 
convince oneself that P' does not vanish if one makes two clusters of m variables. Thus, the minimal 
degree of a polynomial P in 2 to variables that does vanish under C m has degree at least three. In fact, 
it is not too hard to find an expression similar to the one for P' , namely 

P(zi, . . .,Z 2 m) = S[(zi - Zi)(z 2 - Zi){z 3 - Z 4 )\. (54) 

It is not completely obvious that this vanishes under the clustering for to > 2, but one can convince 
oneself that after symmetrization, one indeed does get zero. 

Though it is not going to be easy, one could try to proceed in this way. Constructing the next case, 
namely polynomials that vanish for three clusters of to variables, is already more involved. Writing down 
an explicit form similar to the ones above is not straightforward, but one can for instance symmetrize 
the following combination 


P(zi, . . . , Z 3m ) = S[(zi - Z 5 )(z 2 - Z 5 )(z 3 - z 5 )(z 4 - z 5 ) - {z\ - z 5 )(z 2 - z 5 )(z 2 - z 6 )(z 3 - z 6 )). (55) 


This polynomial is the unique polynomial (up to a constant factor) in 3m variables, of degree and total 
degree 4, that vanishes under formation of three clusters of to variables. We stress, however, that this 
alone does not imply that there are no polynomials of degree three, that vanish under the same clustering 
conditions. 

The lowest degree polynomial for N = 4 and arbitrary to can still be written by symmetrizing an 
expression like the one in Eq. (55), i.e., two terms only, but it seems likely that these expressions become 
more complicated upon increasing N. In addition, having these explicit expressions does not help in 
excluding the existence of lower degree polynomials with the same clustering conditions. 
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VIII. DISCUSSION 

In this paper, we revisited the study of the PES for the v = \/m Laughlin states, in particular the 
rank of the associated reduced density matrix. To determine this rank, we make use of the rank of the 
reduced density matrix for the v = 1 Laughlin state. We showed that to relate the rank for the v =1 /to 
L aughlin state to the case v = 1, one has to prove a bound on the degree of symmetric polynomials 
that vanish under the formation of certain clusters. Though we were not able to finish the proof of 
this statement, we made substantial progress by explicitly constructing a set of polynomials that vanish 
under the clustering, and we proved that the total degree of these polynomials is bounded from below. 

We commented on a possible, though most likely rather hard, route towards finishing the proof. In 
this paper, we concentrated on the Laughlin states. It would be interesting to see if similar methods can 
be used to make progress on different model states, such as the Moore-Read and Read-Rezayi states, 
that exhibit excitations obeying non-Abelian statistics. 


Acknowledgments 

The work of B.M. and E.A. was sponsored in part by the Swedish Research Council. B.M. thanks 
Mohammad Khorrami and B.E. thanks R.Santachiara, A.B.Bernevig, N.Regnault, A.Sterdyniak, P.Zinn- 
Justin, M.Wheeler and D.Betea for discussions. 


Appendix A: Properties of the polynomials r\ 


(*) 


In this Appendix we introduce a family of symmetric polynomials {ri^} defined through the generating 
function 


/ 00 +k \ °° j- 

expf-^pfc-) =J2 r n ] ^ 


\ k =1 / n—0 

The key property of the rk J ’s is their behavior under the clustering transformation C v 


r 

Wi 


(ri x) ) 


= r. 


(m x) 


(Al) 


(A2) 


which is a direct consequence of their definition. Further properties follow from the generating function, 
namely 

• = (-1 ) n n\e n , 

• ri -1) = n\ h n with h n = ^ Xhn m x , 

= —(n — 1 )\p n for n > 1, 

x=0 


• d x Tn ^ 


where the first two relations are obtained by comparison to the generating functions for the e n and h n , 
see for instance 34 . Therefore, this family of polynomials interpolates between power sums p n , elemen¬ 
tary symmetric polynomials e n , and complete homogeneous symmetric polynomials h n . We give one 
additional property, 


r i x +v) 

1 n 



r 


C x)Av ) 

k n—k’> 


(A3) 


that follows from the definition. 
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1. Explicit formulas for r 


(*) 


The generating function can be expanded using Bell’s polynomial^ 39 !, yielding an explicit expression 
for rn \ that is 


r n )= «!53( -iC ) ,(A, ~' (A4) 

\ I 

Ah n 

Alternatively, this expression can be obtained by acting with C x on Newton’s identity expressing elemen¬ 
tary symmetric polynomials in terms of power sums (here, we allow x to be real, and set C x pi = xpi , for 
an infinite number of variables). Another explicit expression is given in terms of a determinant of power 
sums 


r rf } = (-1) 


xpi 1 0 0 

xp 2 xpi 2 ■ ■ : 

: 0 
XPn-l Xp n -2 • ■ ■ Xpi n~ 1 
Xp n Xp n -i ■■■ Xp 2 Xpi 


The first few polynomials are given by 


(A5) 


(x) 

r{ = -xpi, 

(x) 2 2 

*2 = x Pi - X P2, 

r ( 3 x) = —x 3 p\ + 3x 2 p 1 p 2 - 2 xp 3 , 

= x 4 p\ - 6 x 3 plp 2 + 3x 2 pl + 8 x 2 pip 3 - 6 xpi, 

= —x 3 p\ + 10 x 4 plp 2 — I5x 3 pipl ~ 20x 3 p 2 p 3 + 20 x 2 p 2 p 3 + 30x 2 p±p4 — 24a;p5. 


(x) (X) . 

From triangularity it follows that the family {r^ ,... , rA } algebraically spans all symmetric polynomials 
in n variables as long as x 7 ^ 0 . 


2. Degree of r 


(») 


Let us consider the monomial decomposition of rffi 

r n ] =^Z<l { n\ x ) m X- 


(A 6 ) 


Ah r 


We are going to compute the first coefficient, namely q$P(x) for A = (n). If this first coefficient is 
non-zero, the degree of rA J is n. This coefficient is a polynomial of degree n in x , since from (A4) we 
have rn\z, 0 , 0 ,...) ~ (—x) n z n as x goes to infinity. 


We now take x = m to be an integer, and write ri ” 1 ' 1 = C m (r^) = (—1 ) n n\ C m (e n ). By considering 


„(1)N 


the definition of e n , it follows that to > n in order for r 
if the number of variables is less than 1(A). 


(m) 


( z , 0 , 0 ,...) to be non-zero, because m\ = 0 
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Therefore, gi A (a;) vanishes for x = 0,1, • • • , n — 1. It follows from the asymptotic behavior given 
above that 

q (x ^ n)) {x) = (—l) n x(x - l)(x - 2) • • • {x - n + 1). (A7) 

r x \ 

Hence, the degree of r„ isn as long as x 0,1, • • • , n — 1. 


3. Monomial decomposition of r 


(a;) 


In this Appendix we quote the full monomial decomposition of , namely 


(*) = (-!)" 


n 

Ai,..., A r 


n ^ 


(A8) 


i=0 


where A * stands for the transpose of A. The parts of A* are given by A- = /(A) — J2)=i n j( A), or 
equivalently, \\ — A‘ +1 = n,(A). Below, we sketch how this can be established. 


Lemma 1. The monomial expansion of dfr^ 
\i = 0 for i> j + 1. 


a;—0 


involves only partitions A such that A^+i < 0, i.e., 


(x) 

Proof. From the generating function for the rA ; we get 


aW” U = (-*)'"! E III 1 . 

h / 


fcj>l 2 

^1_|-|-fc =77, 


(A9) 


In this expression each term pk i has a monomial decomposition that involve only partitions A with a 
length l( A) < j , from which lemma 1 follows. ■ 


Lemma 2. 77ie monomial expansion of d 3 x r. 
Aj+i < i- 


(x) 


, with i an integer, involves only partitions with 


Proof. The case * = 0 boils down to Lemma [Tj Lemma [5] can be proven by induction on i using 

aj’t’U = " ! E faPL-i-f • 


(A10) 


which follows from taking d J x in (A31, namely 


f)j r {x+y) — ( n ) (y) 

x'n / j \pJ x P n ~P 1 


(All) 


and then choosing x = i — 1 and y = 1. 

Corollary 3. The coefficient qi^\x) is of the form 


x \x ) = - i) A <+i. 


(A12) 
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Proof. Lemma [ 2 ] is equivalent to stating that dfq^\x) = 0 for x = 0,1,..., Aj+i — 1. Thus, x = i is a 
root with degeneracy A- +1 of the polynomial in x q^\x), because A- +1 is the number of parts of A that 
are bigger or equal to i + \ . Since this is true for all i £ N, we have a total of JA \\ = n zeros. Since 
q£\x) is of degree at most n, Corollary [ 3 ] follows. ■ 

Lemma 3. The coefficients c^ are given by 


,W - 


= (-1 y 


n! 


Ai!- - - A„! 


(-l) n 




(A13) 


Proof. The asymptotic behavior of (A6) for x going to infinity yields 


(-1 ) n p? = £ c " A W 

Ah n 


(A14) 


Lemma 3 follows by expanding the left hand side using the multinomial theorem (assuming the number 
of variables p> n) 



E 

k\ ,...,k p >0 
fciH h k p =n 


Tl\ 

h ! - - ■ fcp! 



and then gathering the terms of the r.h.s. into symmetric monomials. 


(A15) 


4. Behavior of rif > under translations 


Translations are well defined in the case of finitely many variables {x±, ■ ■ ■ ,x r }, in which case we set 
Po = r, i.e., the number of variables. In that case L~ is the generator of translations 


l " = E 


d 

dxi 


(A16) 


By Leibniz’s rule, its action on p\ = . p\. is 

L~P\ =^Hj(%-iPA\{j}, Po = r, (A17) 

i>i 

where the nj( A) is the number of parts of A that equal j, and A\ { j } denotes the partition derived from 

(x) 

A by deleting one part that equals j. We can now act on r« , 

L ~r^ = n\ ^7-E ( A18 ) 

\ 1 ^A , 

Ahn j—1 

We can change the summation variable from A to fi = A \ {j}, after noticing that /x is in one-to-one 
mapping with (A,j), since 


j = n-\n\, A = /xU{j}. 


(A19) 
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In particular /(A) = Z(/x) + 1, rij{ A) = rij(fx) + 1 and 3 a = z tl (nj(fi) + l)j. We find 

Mho Z » 

We now split the sum into two parts. First the term |/x| = 77 . — 1 is simply 

, (—(®) 

n\ 2^ -y- PoPn = -xnp 0 r^l v 


Now the reminder is 


/zbn— 1 


1(— , 1 \ (®) 

n - Z _j -“- !>n 1-Im/V = n \ n - l ) r n- 1. 


(A20) 


(A21) 


(A22) 


/xbO 


as can be seen from the determinant formula (A5) by repeatedly developing along the last column until 
the matrix is 1 by 1. These last ‘determinants’ correspond to the factor in the sum above. 

Finally we get 


L r^ x) =n(n- 1- xp 0 )r ( n x 2 1 


(A23) 


from which Eq. (52) in the main text follows. 


Appendix B: An alternate set of generators 


In this Appendix, we briefly introduce an alternate set of generators p[ , that could be used in the proof 
in section [v] This set of generators of A q=m N is constructed to satisfy p' = pi for * < IV, and C m (p') = 0 
for i > N, and have total degree i. In the construction, to keep track of the number of variables of the 
power sums, we introduce an additional index IV, so 


Pi,N(xi, ... ,x N ) = x\ H- \-x l N . (Bl) 

Recall that An = {pi,n, • • • ,Pn,n} are algebraically independent and generates all symmetric poly¬ 
nomials in N variables. In particular for each i € IN* there exists a unique polynomial Ti imj Ar in N 
variables, such that 

mp it N = T iim ,iv(rapi,jv, • • • at). (B2) 


The p' are now defined as follows 


Pi,q = 


Pi,q , 1 < i < N 


• ■ ■ iPN^q)) N <C 2 ^ Q. 


(B3) 


By construction, they obey C m (p' ) = 0 for i > N+ 1, are non-zero, and have total degree i. In addition, 
they form an alternate generating set A' q = {p'i t q, ■ ■ ■ tP' q ,q\ °f A q , because the pi tq of the generating set 
A q can be expressed in terms of the p' i q in Eq. ( |B3[ ). 

Since, as is shown in Corollary [lj pn+ 1 is the unique (up to a scale factor) polynomial that vanishes 
under the clustering C m , we find that p' N+1 oc Pn+i, and that L~p' N+1 = 0. In addition, by a direct 
calculation, one finds that L~p' i q = ip[_ l q for 1 < i < N and N + 1 < i < q. 

Namely, by acting on both sides of the definition of Ti^n, Eq. (B2) with L~ gives 
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N 

L~T itrrit N(mpi t N, ■ ■ ■, mpN,N) = qTi, m ,N-,\{mpi,N, ■ ■ ■ ,nipN,N) + tnjPj-i,NTi tm ,N.j(mpi t N, ■ ■ ■ ,mpN,N ) 

j =2 

= impi-i >N = iTi- 1:myN (mpi : N, ■ ■ ■ ,mpN,N), 

where T, , m ,N;j denotes the derivative of Tj )TO> jv with respect to its jtli argument. In particular by setting 
Xj = mpj t N, we find that for general arguments, 


N 

iTi-i } m,N{Xi,... ,Xn) = qTi im ^-\{X\,... ,Xn) + ^2,jXj-\Ti, m .,N;j{.Xi,... ,Xn)- (B4) 

j=2 


We can now act with L on both sides of the definition of p' , 
we find that, for i > N + 1, 


Eq. (B3|. Using the relation Eq. (B4), 


L Pi,q= L Pi,q~ L Ti(p liq ,...,p N , q ) 

= iPi-l >q - . . . ,p N ,q) (B5) 

= iPi-l,qi 


which is what we wanted to show. We note that in the case that i = N + 1, the only thing that changes 
in the argument above is that the left hand side of Eq. (B4) is replaced by iX n, leading to the result 
L~p ' N+ 1 = 0, which we showed in the main text using a different method. 

Finally, we mention that it is also possible to prove that the degree of p' equals i for A^+l < i < 2AT+1, 


directly from the definition. We believe that the degree of p' i q also equals i for 2N + 1 < i < q, but did 
not find a proof of this. 
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